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The relation of the recently proposed E(5) critical point symmetry with the 
interacting boson model is investigated. The large-N limit of the interacting 
boson model at the critical point in the transition from U(5) to 0(6) is ob- 
tained by solving the Richardson equations. It is shown explicitly that this 
algebraic calculation leads to the same results as the solution of the Bohr 

X 

■ differential equation with a potential 
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The study of phase transitions is one of the most exciting topics in Physics. Recently the 
concept of critical point symmetry has been proposed by Iachello [1]. These kind of sym- 
metries apply when a quantal system undergoes transitions between traditional dynamical 
symmetries. In Ref. [1] the particular case of the Bohr Hamiltonian [2] in Nuclear Physics 
was worked out. In this case, in the situation in which the potential energy surface in the 
ft-'y plane is 7-independent and the dependence in the ft degree of freedom can be modeled 
by an infinite square well, the so called E(5) symmetry appears. This situation is expected 
to be realized in actual nuclei when they undergo a transition from spherical to 7-unstable 
deformed shapes. The E(5) symmetry is obtained within the formalism based on the Bohr 
hamiltonian, but it has also been used in connection with the Interacting Boson Model 
(IBM) [3]. Although this is not the form it was originally proposed [1], it has been in fact 
argued that moving from the spherical to the 7-unstable deformed case within the IBM one 
should reobtain, at the critical point in the transition, the predictions of the E(5) symmetry. 
This correspondence is supposed to be valid in the limit of large number N of bosons, but 
the calculations with the IBM should provide predictions for finite N as stated in Ref. [4]. 
In this letter, on one hand we calculate exactly the large N limit of the IBM at the critical 
point in the transition from U(5) (spherical case) to 0(6) (deformed 7-unstable case). On 
the other hand, we solve the Bohr differential equation for a ft 4 potential. Both calculations 
lead to the same results and are not close to those obtained by solving the Bohr equation for 
an infinite square well (E(5) symmetry). We also show with two schematic examples that 
the corrections arising from the finite number of bosons are important. With this in mind, 
the IBM calculations still provide a tool for including corrections due to the finite number 
of bosons. 

In Ref. [1] the Bohr Hamiltonian is considered for the case of a 7 independent potential, 
described by an infinite square well in the ft variable. In that case, the hamiltonian is 
separable in both variables and if we set 

*{p,iA) = f{PMiA) (i) 
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where 9i stands for the three Euler angles, the Schrodinger equation can be split in two 
equations. The solutions of the (7, 9j) part were studied in Ref. [5] and tabulated in Ref. [6]. 
Iachello solved the (5 part and found that the f(/3) functions are related to Bessel functions. 
The main results are illustrated in Table I and Fig. 1 of Ref. [1] . These results are obtained 
from a geometrical picture and we would like to investigate its relation with the interacting 
boson model. 

The geometrical interpretation of the abstract IBM hamiltonian can be obtained by 
introducing a coherent state [7-9] which allows to associate to it a geometrical shape in 
terms of the deformation variables ((3, 7) . The basic idea of this formalism is to consider 
that the pure quadrupole states are globally described by a boson condensate of the form 

\g;N,p, 1 ) = -±=(Tl) N \0) , (2) 

where the basic boson is given by 

(3) 
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+ (3 cos 7c4 H — t=/3 sin j(dl + dL 2 ) 
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which depends on the (3 and 7 shape variables. The energy surface is defined as 

E N (P n ) = {g;N,P n \H\g;N,P n ) , (4) 

where H is the IBM hamiltonian. Here we are interested in the case in which the hamiltonian 
undergoes a transition from U(5) to 0(6) and, consequently, the corresponding potential 
energy surfaces are 7-independent. 

In order to investigate the geometrical limit of the IBM in the transitional class going 
from U(5) (spherical) to 0(6) (deformed 7-unstable) the most general (up to two-body 
terms) IBM hamiltonian is, 

H = e d h d + K P^P + «iL • L + k 2 Q x=0 ■ Q x=0 + k 3 T 3 ■ f 3 + k a T a ■ f 4 (5) 

where fi d is the d boson number operator, and 
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pt = I (dt.dt_ s t. s t) > 



i = viois x j)«, 

gx=o = ( s t X d + Sx 5) (2 \ 
f 3 = (Jxdf\ 
f 4 = (dt x J)W 



(6) 
(7) 
(8) 
(9) 
(10) 



The scalar product is defined as T L -f L = J2m(~^) M ^lmTl-m, where T LM corresponds to the 
M component of the operator Tl- The operators d m = (— l) m d_ m and s = s are introduced 
to ensure the correct tensorial character under spatial rotations. The corresponding energy 
surface is obtained from Eq. (4) 



E(N,P) 



N 



+ 



1 + /3 2 
N(N - 1) 
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5k 2 + (3 2 (e d + 6/ti + k 2 + -k 3 + -K4 



(1 + /9 2 ) 2 

The condition to find the critical point is 



(l-(3 2 ) 2 Ao2 18 oA 
4 oo 



(11) 



(d 2 S(JV, ) 9)/d^) /jM) = 



(12) 



and gives the following relation among the hamiltonian parameters 
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e d = -6/Ci + 4k 2 - -/c 3 - -/c 4 + (A/" - l)(/c - 4/c 2 ). 

O 



(13) 



Thus the most general energy surface at the critical point in the U(5)-0(6) phase transition 
is 



E crit (N, (3) = 5Nk 2 + N(N - 1) 



K ° , ( a , 18 ^ ^ 



(14) 



These expressions are consistent with those obtained in Ref. [10] for a slightly different 
hamiltonian. Note that (14) completely defines the form of the potential up to a scale and 
an energy translation. The expansion of this critical energy surface around (5 = is 



E crit (N, (3) « 5k 2 N + 'jNiN - 1) + N(N - 1) (kq - Ak 2 + ^k 4 ) * 1 - 2 >" ; 



(15) 



whose leading term is f3 4 . Alternatively, one can carry out the transformation (3 2 /(l + {3 2 ) — > 
(3 2 and finds f3 A as the critical potential. 

In order to make some calculations to illustrate the large N limit in the IBM at the 
critical point in the U(5)-0(6) phase transition and the corresponding finite N corrections, 
we propose two schematic transitional hamiltonians. The first one is 

^ = xfld + WTi p]p ■ (16) 

The corresponding energy surface is obtained from Eq. (11) with e d = x, n = and all 
the rest of the parameters equal to 0, 



Ej(N,P)=N 



(3 2 1-x fl-(3 2 \ 2 

X T77 + 



(17) 



1+/9 2 4 \1 + (3 2 J 

The condition to localize the critical point, Eq. (13), gives in this case x\ = 0.5. In Fig. 1 
we represent as an example the energy surfaces for the hamiltonian (16) (left panel) with 
three selections for the order parameter x: one at the critical point, one above that value 
and one below it. For x > x c an equilibrium spherical shape is obtained, while for x < x c 
the equilibrium shape is deformed. The value x c gives a flat (3 4 surface close to (3 — 0. 
The second schematic hamiltonian we propose is 

Hn = xh d - ^Q x=0 ■ Q x=0 , (18) 

The corresponding energy surface is obtained from Eq. (11) with e d — x, k 2 — —^jf- and 
all the rest of the parameters equal to 0, 

En(N, P) = -^ + P 2 )l^ + Nx Y ^ W2 - A(N - 1)(1 - x) . (19) 

Condition (13) gives in this case the critical point x 1 / = f^Ef that in the large N limit gives 
4/5. 

In Fig. 1 the corresponding energy surfaces are plotted in the right panel. Same com- 
ments as in the preceding case are in order. Thus, we conclude that, in the transition from 
spherical systems to 7-unstable deformed ones, the critical point in IBM should be associated 



to a j3 4 potential rather that to an infinite square well. The question is then how different 
are the E(5) predictions from those obtained with a f3 A potential? In order to investigate 
this point we have solved numerically the Bohr hamiltonian for a potential (3 A . The results 
for energies are presented in Table I and in Fig. 2. Here we keep the label £ used in the 
E(5) case. It is related to the label rip = sometimes used in the U(5) classification, 

by np — £ — 1, where is the U(5) label and r is the 0(5) label. Particularly interesting 
are the energy ratios given in Table II which have been used in recent works to identify 
possible nuclei as critical. In this table the E(5) and (3 4 values are shown for comparison. 
The labeling of the states is L^ T . 

Besides the excitation energies, B(E2) transition probabilities can be calculated using 
the quadrupole operator 



If 2 ' = t (3 pgftfc) cos 7 + ^ (^M) + Z><2 2 (0O) sin 7 



(20) 



where t is a scale factor. In Table II two important B(E2) ratios are given for E(5) and (3 4 
cases. In Fig. 2 the B(E2) values for a f3 A potential are shown besides the arrows. They are 
given normalized to the B(E2; 2^ — > 0i )0 ) value which is taken as 100. 

Comparing Figs. 1 and Table I in Ref. [1] with the present Fig. 2 and Table I we 
can observe important differences between E(5) and (3 A potentials. In order to see which 
is the actual large N limit of IBM we have performed calculations with the IBM codes for 
hamiltonians Hi (Eq. 16) and H n (Eq. 18) at the critical point for different number of 
bosons. These codes allow to manage a small number of bosons, typically 20. In Fig. 3 the 
results of these calculations are shown with a full line for Eq. (16) and with a dashed line 
for Eq. (18). The values for E(5) and (3 A potentials are shown as dotted lines as references. 
The last two panels labeled with Ri and i?2 refer to the B(E2) ratios presented in Table II. 

From Fig. 3 it is clear that the finite N effects are important and depend on the precise 
form of the hamiltonian used. However, it is difficult to conclude whether E(5) or (3 4 is the 
large N limit of the corresponding IBM hamiltonian. It is necessary to perform calculations 
with larger values of N. Fortunately, Dukelsky et al. [11] have recovered an exactly solvable 



model for pairing proposed by Richardson in the 60's [12]. Following Ref. [11] we have 
solved the Richardson's equations and obtained the exact eigenvalues for the hamiltonians 
(16) and (18) up to N = 1000, so approaching the large N limit of the corresponding IBM 
hamiltonians. Details of this method will be given in a longer publication. In Fig. 4 we 
present the results of these calculations for energy ratios up to iV = 1000 and B(E2) ratios 
up to N = 40 together with the corresponding values for the E(5) symmetry and the f3 4 
potential. From this figure it clearly emerges that the large N limit for the studied IBM 
hamiltonians corresponds to the (3 4 potential. Both hamiltonians Eq. (16) and Eq. (18) 
converge to the same results in the large N limit, although the corresponding corrections for 
finite N are quite different (see Fig. 3). 

We conclude that the large N limit of the IBM hamiltonian at the critical point in 
the transition from U(5) (spherical) to 0(6) (deformed 7-unstable) is represented in the 
geometrical model by a (3 4 potential. The results are similar but not close to those of an 
infinite square well as in the E(5) critical point symmetry. The analysis of the IBM energy 
surface followed by an IBM calculation, as presented in Ref. [13], can provide the appropriate 
finite N corrections and thus lead to the identification of nuclei at the critical points. In that 
work a systematic study of the properties of the Ru isotopes allowed to select the appropriate 
form of the hamiltonian. Once it is fixed the construction of the energy surfaces identify 
the critical nucleus ( 104 Ru in that case). The corresponding IBM calculation for the critical 
nucleus then provides the correct finite N corrections. We believe that this is a fundamental 
step if we wish to robustly identify the spectroscopic properties that signal the presence of 
criticality in the atomic nucleus. 
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TABLES 



TABLE I. Excitation energies for a [3 A potential relative to the energy of the first excited state. 





£ = i 


e = 2 


£ = 3 


e = 4 


T = 


0.00 


2.39 


5.15 


8.20 


T = 1 


1.00 


3.63 


6.56 


9.75 


T = 2 


2.09 


4.92 


8.01 


11.34 


T = 3 


3.27 


6.26 


9.50 


12.95 



TABLE II. Energy and B(E2) transition rate ratios in the E(5) symmetry and for the /3 4 
potential. 





^41,2/^21,1 E 02 /E 2lA E 0l3 /E 2lA E 02 /E Ql3 


p B(E2;4i, 2 ^2i,i) p B(£2;0 2 ,o->2i,i) 
^1 - B(E2;2i,i-^0i, ) ^ ~ B(£2;2i,i-Oi,o) 


E(5) 

/3 4 


2.20 3.03 3.59 0.84 
2.09 2.39 3.27 0.73 


1.68 0.86 
1.82 1.41 
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FIGURES 




FIG. 1. Representation of the energy surfaces for N = 20 as functions of the shape parameter 
(5 obtained for two schematic hamiltonians, Eq. (16) (left panel) and Eq. (18) (right panel). In 
each case three values of the order parameter are presented, one at the critical value, one above and 
one below that value. The curves have been arbitrarily displaced in energy so as to show clearly 
the behavior. 
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FIG. 2. Schematic spectrum for a (5 potential. Numbers close to the arrows are B(E2) values. 
These are relative to the transition 2i ; i — > O^o whose B(E2) value is taken as 100. 
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5 10 15 20 
N 

FIG. 3. Variation with the number of bosons (up to N = 20) of selected energy and B(E2) 
ratios for IBM calculations performed at the critical points of hamiltonian (16) (full line) and (18) 
(dashed line). The corresponding E(5) and (3 values are marked with dotted lines. 
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FIG. 4. Same as Fig 3 but here the number of bosons runs up to 1000 in the energy ratios and 
up to 40 in the B(E2) ratios. 



13 



